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Abstract 



We study linearized gravitational perturbations of extreme black hole solutions of 
the vacuum Einstein equation in any number of dimensions. We find that the equations 
qh governing such perturbations can be decoupled at the future event horizon. Using these 

equations, we show that transverse derivatives of certain gauge invariant quantities blow 
bD up at late time along the horizon if the black hole solution satisfies certain conditions. We 

find that these conditions are indeed satisfied by many extreme Myers- Perry solutions, 
including all such solutions in five dimensions. 

>; 

CO 

ov 1 Introduction 

Extreme black holes have theoretical importance in understanding of quantum theory 
of gravity. For example, Bekenstein-Hawking entropy of supersymmetric black holes was 
explained by counting BPS states in the view of string theory [T]. Furthermore, a duality 
called the Kerr/CFT correspondence between extreme black holes and a two-dimensional 
conformal field theory was proposed |2], and the entropy of the black holes was reproduced 
5-j | as the statistical entropy of the dual CFT. 

Recently, it was shown that extreme Reissner-Nordstrom and Kerr black holes are classi- 
cally unstable against test scalar field perturbations [3H6]. Subsequently, the proof is extended 
to all other extreme black holes [TJ. They showed that the second transverse derivative blows 
up at the horizons as d^cf) ~ v, where <fi is the scalar field and we took the ingoing Eddington- 
Finkelstein coordinates, (v,r). For extreme Kerr black holes, the similar instabilities were 
also found in gravitational and electromagnetic perturbations [TJ. 

We arise following questions: Are all extreme black holes unstable against gravitational or 
electromagnetic perturbations? If not, what is the condition for the instability? In this paper, 
we address these questions by studying the perturbations of any extreme black holes. We use 
Geroch-Held-Penrose (GHP) formalism in higher dimensions developed in Refs. j8,9j to study 
the perturbations. So, in sectionj2l we give a brief review of the GHP formalism. We introduce 
gravitational, electromagnetic and scalar field perturbation equations based on the formalism. 
They are regarded as higher dimensional analogues of Teukolsky equations although they are 
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not decoupled equations for gravitational and electromagnetic perturbations in general. In 
section^ we introduce the most general expression for extreme black hole and express them 
in the view of the GHP formalism. In sectionJU we study the scalar field perturbation. 
Although the scalar field perturbation on any extreme black holes has been already studied 
in Ref. [7], we revisit the problem using the GHP formalism. We find that all extreme 
black holes are unstable against scalar field perturbations as shown in Ref. [7J. In sectionJU 
we study electromagnetic perturbations. We find that, near the horizon, electromagnetic 
perturbations satisfy decoupled equations. Using the decoupled equations, we show that the 
perturbations do not decay along the future event horizon if a certain operator on the horizon 
has a zero eigenvalue. In sectionJHl we study gravitational perturbations. By the similar way as 
electromagnetic perturbations, we can show the non-decay of the gravitational perturbations 
if a horizon operator has a zero eigenvalue. In addition to that, if the background geometry 
is algebraically special, the first or second transverse derivatives of the perturbation variables 
blow up along the horizon. The eigenvalues for the horizon operators have been calculated for 
some extreme black holes. In sectionJTJ we see that there are zero eigenvalues in the horizon 
operators for all higher dimensional extreme black holes with zero cosmological constant as 
far as we calculated. The final section is devoted to discussions. 



2 Geroch-Held-Penrose formalism in higher dimensions 

We study the perturbation of the general extreme black holes using the Geroch-Held- 
Penrose (GHP) formalism in higher dimensions developed in Refs. [8j[9]. In this section, we 
give a brief review of the GHP formalism. In the formalism, we use a null basis {eo, e±, e{\ = 
{£, n, rrii} (i = 2, • • • , d — 1) which satisfies 

£ 2 = n 2 = £ ■ rrii = n ■ rrii = , £ ■ n — 1 , m ; • rrij = Sij . (2-1) 

We define the covariant derivatives of basis vectors as 

L ab = V 6 4 , N ab = V b n a , Mi b = V b m ia , (2.2) 

and 

Pij — , Ti = Ln , Ki = Lio . (2-3) 

The orthogonal relations (12. ip are invariant under spins, boosts and null rotations defined as 
follows. Spins are local SO(d — 2) rotations of the spacial basis {rrii}: 

rrii -> Xijrrij , (2.4) 

where G SO(d — 2) depends on the spacetime coordinate x^. Boosts are local rescaling 
of the null basis: 

£ -> X£ , n ->• n/X , (2.5) 
where A is any real scalar function. Null rotations about I and n are 

I — y £, n — > n + Zifrii — z 2 £/2 , m ; — > m ; — Zil , (2.6) 
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and 

+ z'frii — z /2 n/2, n — > n , m 8 — >■ rrii — z^n , (2.7) 

where Zj and z[ are real functions of x^. 

In the GHP formalism, we maintain the covariance with respective to spin and boost 
transformations. An object T^...^ is a GHP scalar of spin s and boost weight 6 if it transforms 
by the spins and boosts as T^...^ — > Xi 1 j 1 ■ ■ ■ Xi s j s Tj 1 ...j s and T^...^ — » A^T^...^. For example, 
the quantities p^, Tj, are GHP scalars with 6 = 1, 0, 2, respectively. We also define priming 
operation: T^...^ — > t , where T( ^ is the object obtained by exchanging £ and n in the 
definition of T^...^. 

We define GHP scalars obtained from Weyl tensor C a bcd as 

&ij = CoiOj j — Cuij , (2-8) 

^ijk = Coijk , ^ijk = Cujk , = Coioi , = Cioii , (2.9) 
% = , ®m = C m , $ = Coioi , $g = %) , = %] , (2.10) 

where f2, $, and f2' have boost weights b = 2,1,0,-1,-2, respectively. The null 
vector i is called multiple WAND (Weyl-aligned null direction) iff all boost weight +2 and 
+ 1 components of the Weyl tensor vanish. The spacetime admitting the multiple WAND is 
called algebraically special spacetime. We can also obtain GHP scalars from Maxwell field 
strength F a b are 

<p i = F oi , F — Foi , Fij = Fa , p\ = F u , (2.11) 

where tp, F and tp' have boost weights b = 1, 0, —1, respectively. 

The partial derivatives of GHP scalars, such as i^d^T^...^, n^d fJL T il ...i s or m^c^T^...^, are 
not GHP scalars. It is convenient to define derivative operators which are covariant under 
spins and boosts as 

s 

!>'/;>..,. = ^T h ... is - bL 10 T h .. 4s + ^Mf r0 T n ... lr _ lfclr+1 ... Js , (2.12) 

r=l 

s 

\>' T h-i s = n»d,T n ... is - 6L n T n .., s +J2 M ti T ^r- l ki r+ v--i s , (2-13) 

r=l 
s 

^T n ... js = m»d,T n .., s - bL u T n ... 3a + j- • (2.14) 

r=l 

They are called GHP derivatives. We can check that bT^...^, b'T^...^ and <5iTj 1 ...j s are all GHP 
scalars, with boost weight (6 + 1,6 — 1,6) and spins (s, s, s + 1). 

The GHP scalars defined above are not independent because of Ricci equations, [V M , V„] V p = 
RftvpaV 7 , Bianchi equations, V \\C — 0, and Maxwell equations, dF = d* F = 0. The re- 
lation for the GHP scalars in Einstein spacetimes = Ag^ are summarized in appendixICl 
Since these equations are invariant under the spins and boosts, they are written by GHP 
scalars and their GHP derivatives. 

In the GHP formalism, the Klein-Gordon equation (V 2 — p 2 )0 = is written as 

(2b'b + SA + p'p - 2rA + pp' - p 2 )</> = . (2.15) 
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From appropriate linear combinations of equations in appendixJCj we can obtain useful equa- 
tions for studying electromagnetic and gravitational perturbations [9]. They are written as 

2d — 3 

(2p'p + 3,3, + p'p - Arfo + $ - -j—[ A ^ + (~ 2t ^ + 2r A + 2$ S + 
= [«!>' + p3 + (3p) + Q>'k) + pr + K p' + tt] F 16 ^ 
+ (pb' + kk' + pp')(/> + («8 + p 2 + kt + , 

and 

(2b'b + 3 fc 3 fc + pb - 6r fc 3 fc + 4$ - -— yA)% 

+ 4(r fc 3 (i - r (i 3 fc + $g, fe + 4$^)%) + 2$ ifci jfi w + 4«; fc b'(^ fo> + 
= [p3 + rp + r'p + up' + Q)'k) + (3p) + + p 2 $ + Kp^' + (pb'fi + pp' + . 

(2.17) 

The right hand sides of these equations are very long so we wrote them schematically. In this 
paper, we do not need the detailed expressions of the right hand sides. 



3 Extreme black holes in the GHP formalism 

We consider general extreme black holes. The metric of the extreme black holes can be 
written as [10] 

ds 2 = L 2 (x)[—r 2 F(r, x)dv 2 + 2dvdr] + 7 Q/ 3(r, x)(dx a — rh a (r, x)dv)(dx^ — rh?{r, x)dv) , (3.1) 

where functions (F,j a p,h a ) and L 2 are smooth function of {r, x a } and {x a }, respectively. 
The horizon of the spacetime is located at r = 0. In the metric, there is a residual coordinate 
transformation, r — > T(x)r. We choose the free function T(x) so that F(r = 0, x) — 1 is 
satisfied. In this paper, we focus only on Einstein spacetimes satisfying R^ v = Ag^. 

We assume that the background metric have n rotational symmetry generated by d/dcp 1 
(I = 1, 2, • • • , n). Then, the metric can be written as 

ds 2 = L 2 (y)[-r 2 F(r,y)dv 2 + 2dvdr] + 'j AB (r,y)(dy A - rh A (r,y)dv)(dy B - rh B (r,y)dv) 
+ 2j AI (r,y)(dy A - rh A (r,y)dv)(d(p I - rh I (r,y)dv) 

+ 1ij(r, y)(d(f) 1 - r//(r, y)dv)(d(j) J - rh J (r, y)dv) . (3.2) 

We impose further assumption on metric functions as 

h A (r, y) = 0(r) , h z (r, y) = k 1 + 0{r) , lAI (r, y) = 0(r) , (3.3) 

where k 1 are constants. These assumptions are true for a large class of extreme black holes [TlT 
[T5] . Under these assumptions, the near horizon geometry of the metric (13.21) takes "standard" 
form: 

ds 2 = L 2 (y)[-R 2 dV 2 + 2dVdR]+-f AB (y)dy A dy B + -f IJ (y)(d(j) 1 -rk I dv)(d(j) J -rk J dv) , (3.4) 
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where we took the double scaling limit: r = eR, v — V/e and e — > 0. The induced metric on 
the horizon is written as 



ds 2 H = g^dsfdx" = iA B (y)dy A dy B + llJ (y)d<p I d<p J . (3.5) 
We take null basis {eo,ei,ej} = {£,n,rrii} in the general extreme black hole metric ( 13.11) 

as 

2 r 2 F 2rh i 1 
£ = yd v + —j—d r + — —it , n = —d r , m i = e i , (3.6) 

where h % = /i Q e„ and e* is an appropriate orthogonal basis for 7 aj g. The null basis (13. 6p is 
regular at the future horizon r = 0. Using the basis, we can obtain GHP variables. The full 
expression of the GHP variables are summarized in appendixfX] Here, we focus on and Ki 
since they will be important later. They are given as 

2r 2 F ■ 2r 2r r 2 F 

^ = L '* , Pij = -j-h a i(j e i)a + — /i fc ey(ei )Q ) >fc + — ^(e^)' , (3.7) 

where ^ = ef<9 M . We chose the residual gauge freedom so that F(r — 0,x) — 1 is satisfied. 
Thus, we obtain Fj = 0(r). Therefore, we have Ki = 0(r 3 ). From the assumption (13.31) . 
h a \ r= o is constant and we have h a j = 0(r). Thus, the first term in pij is 0(r 2 ). In the second 
term, there is a derivative operator, h k dk = ^d^i + 0(r). Since the d/d^i is a Killing vector, 
its operation to background variables vanishes and the second term is also 0(r 2 ). The last 
term is trivially 0(r 2 ). Therefore, we can conclude that is second order in r. By the 
similar way, we obtain the near horizon expression of the GHP variables as 

Pi 3 = 0{r 2 ) , «, = 0(r 3 ) , r, = h + 0(r) , 

4, = 0(1), «J = 0, 7* = 0(1), 
L io = T + 0(r 2 ) , L n = , L u = ^ + O(r) , 
Mj = (9(r 2 ) , M) x = 0(1) , M) k = 0(1) . 

Components of Weyl tensors with boost parameter +2 and +1 are given by Newman- Penrose 
equations flNPTj) and flNP2j) as 

% = -tPij + 3j K i - PifcPfcj - k^t] - TiKj = C(r 3 ) , 
^ijfc = 2(r l p [ifc] + K^p'p] - 8 b -|Pi|jfe]) = C(r 2 ) . 

Now, we consider Components of Weyl tensor with boost parameter 0. In Ref. [8], it was 
shown that the induced Riemann tensor on a spacelike surface which is orthogonal to null 
vectors I and n is written as 

2A 

R W = 2 Pk[i\p'i\j] + 2p'k{i\Pi\j] + ®ijki + ^7— [%|fe% • ( 3 - 10 ) 
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Thus, we have 

^ijkl Hijkl 



^ (3 11) 

*g = ~{Ra ~ J^*8*) + 0(r) , $ = -\{R - (d ~^~ 3) A) + O(r) , 

where -Rjjfcz is the induced Riemann tensor on the horizon (j3.5p . From the antisymmetric part 
of Eq. (INFi|) . we obtain 



1 ;dk + — Ld£ 2 A k 



4L 2 4L 2 



+ 0(r). (3.12) 



The GHP derivatives are 



\>T h ... is = j{d v + r^df - b)}T h ... is + 0(r 2 ) 
1 s 

P ^ii— i a — 7TjdrTi 1 ...i s + ^ ^ -^i r l^ii— if-ifcir+l—t 



(3.13) 



r=l 



where V is a covariant derivative with respect to the horizon induced metric (13.51) . 

4 Scalar field perturbations 

4.1 Conserved quantity on the horizon 

First, we consider the scalar field perturbation equation ( 12. 15ft . The instability of the 
massless scalar field perturbation on any extreme black holes has been already shown in Ref. 
[7]. Here, we revisit the problem including the massive scalar field using the GHP formalism. 

Substituting near horizon expressions of GHP variables and derivatives 03.8H3.l3]) into the 
Klein-Gordon equation (12. 151) . we have the scalar field equation near the horizon as 

8 V [2L(2b'0 + P '<f>)] = A <j> + 0(r) . (4.1) 

where the operator Ao is defined ad3 

A 4> = -Vi(L 2 Vi4>) + itfmri + fi 2 L 2 (p , (4.2) 

where we decomposed {</> 7 }-dependence of <p by Fourier modes e imi ^\ that is, di<f) = mj<j). 
We will do same decompositions for electromagnetic and gravitational perturbations. The 



1 For axisymmetric perturbations, this operator Ao coincides with the operator defined in the study of 
perturbation of near horizon geometries |16j . We will also define operators A s for electromagnetic (s = 1) and 
gravitational (s = 2) perturbations by similar ways as the scalar field, which relate to defined in Ref. [TC] 
as C?W = Ai and = A2 + 2 for axisymmetric perturbations. 
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derivation of the equation is given in appendixjBj Hereafter, we focus on axisymmetric 
perturbations mj = 0. For axisymmetric perturbations, the operator Ao is self-adjoint, 
(Y 1 ,AoY 2 ) = (-Ao^i, ^2)) with respect to an inner product 

(Y U Y 2 )= f Y*Y 2 , (4.3) 

J H 

where f H = f d d ~ 2 xy/g and g is the determinant of the horizon induced metric defined in 
Eq.tfEi]). 

We assume that operator Ao has a zero eigenvalue. This is always true for massless case 
fi = since we have AqY = when Y is a constant. For massive scalar field, Ao can also have 
zero eigenvalues depending on the value of /1 2 and background geometries. We will discuss 
the existence of the zero eigenvalues in sectionJ772J We denote the eigenfunction for the zero 
eigenvalue as Y. Operating (Y, *) to both side of Eq. fl4.ip . we obtain 

^ = , Jo = Jy* [2L(2b> + p'<j>)] . (4.4) 

where we used (Y, Ao4>) = (AqY, <fi) = 0. Hence, Iq is a conserved quantity along the horizon. 
Note that the GHP derivative b' contains only the radial derivative d r . (See Eq. ( 13.131) .) Thus, 
the integrand in the Jo is written by the linear combination of <9 r and (p. Therefore, we can 
conclude that if Jo 7^ at an initial surface, d r <p and do not both decay along the future 
horizon as v — > 00. 

4.2 Instability against scalar field perturbations 

We assume that (f> and its angular derivatives d/dy A decay along the horizon. For extreme 
Reissner-Nordstrom and Kerr black holes, it was shown that <j) decays along the horizon. So, 
this assumption seems likely also for any other extreme black holes. Then, at the late time, 
the conserved quantity J approaches 

J ~ 2 ! Y*d r (p , (v 00) . (4.5) 
Jh 

Now, we differentiate Eq. (l2.15p by r. Near the horizon, the equation can be written as 

d v [2L«9 r (2b'0 + p»] = (A - 2)d r <p + V 4> + 0{r) , (4.6) 

where, in the linear operator T> , there is no radial derivative d r . Thus, we have T> <f) — > as 
v — > 00. The derivation of above equation is in appendixjBj Operating (Y, *) to both side of 
Eq. ()4.6p . we obtain 

^ ~ -Jo , (v ^ 00) , Jo = f Y* [2L«9 r (2b'0 + p'<f>)} (4.7) 

Therefore, the quantity J blows up linearly in time v as 

J ~ -I v , (v -> 00) . (4.8) 

The integrand in the Jo is written by the linear combination of d^(f), d r <f) and <fi. Since we 
assumed that <fi decays at the horizon, either d^cf) or d r cf) blows up at the horizon. This implies 
the instability of the extreme black holes against scalar field perturbations. 
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5 Electromagnetic perturbations 



5.1 Decoupled equation on the horizon 

Secondly, we consider electromagnetic perturbations. We consider the Maxwell field as 
test field and, thus, it vanishes in the background. It follows that GHP scalars obtained from 
Maxwell field strength are all invariant under the infinitesimal coordinate transformations and 
basis transformations ( I2.4j) - (12.7j) . The number of physical degrees of freedom of the Maxwell 
field in ^-dimensions is d — 2. The number of components of tfi is also d — 2. Thus, we can 
expect that tpi has all physical degrees of freedom of electromagnetic perturbations and it 
would be nice if we can obtain decouple equations for tpi. The right hand side of Eq. (12. 16j) 
contains coupling terms between tpi and other components of the perturbation. So, ifi does 
not decouple in general. We can see that, however, all the terms in right hand side of Eq. (j2. 16j) 
are multiplied by p, k, Q or \F From Eqs. (l3.8p and (13. 9p . they are at most 0(r 2 ). Thus, the 
right hand side is 0(r 2 ). Therefore, in Eq. (l2.16p and its radial derivative, the right hand side 
is zero at the horizon. 



5.2 Conserved quantity on the horizon 

We consider left hand side of Eq. (12. 16j) neglecting the right hand side. Using near horizon 
expressions of GHP variables and derivatives (I3.8H3.13]) . we obtain 

d v [2L(2V<Pi + P'^i)] = Am + 0{r) . (5.1) 
where we define the operator A± as 

Ai<Pi = -—V j{L A V j(pi) + ( 2 + 3zfc 7 mr - TT^^ J <Pi + L * + \^9ij) <fj 

1 . 1 \ (5 ' 2) 

-{dk)ij + 2(k - diL 2 ))^ - —{dL%kA ifj . 

The derivation of the equation is written in appendixjB] Hereafter, we focus on the axisym- 
metric perturbations mi = 0. We define an inner product as 

(Y\Y 2 )= [ L\YtyY? . (5.3) 

J H 

For axisymmetric perturbations, we have (F 1 , AlF 2 ) = (AiY 1 ^ 2 ), that is, the operator A\ 
is self-adjoint. We assume that the operator A\ has zero eigenvalue. Although the existence 
of the zero eigenvalue is not obvious, we will see that many extreme black holes satisfy this 
assumption in section ]T72l We denote the eigenfunction by Yj. Operating (Y,*) to Eq. (15 . 1 1) . 
we have 

^ = o , h = j y;[2L 3 (2 P v. + Ml , (5.4) 

where we used (Y,A\ip) = (AiY,ip) = 0. Therefore, I\ is a conserved quantity along the 
horizon. Thus, If Ii ^ at an initial surface, d r (fi and ipi do not both decay along the future 
horizon. 
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We assume that tpi and its tangential derivatives along the horizon decay as v — > oo. Now, 
we differentiate Eq. fl2.16p by r. Then, the right hand side becomes 0(r) and still vanishes on 
the horizon. Thus, we have 

d v [2Ld r (2V<Pi + (/<fi)] = + V m + 0(r) , (5.5) 

where, in the linear operator T>i, there is no radial derivative d r . Hence, we have 1>\^Pi —t 
{v — > oo). The derivation of the equation is in appendix|B] Operating (Y, *) to above equation 
and taking limit of v — > oo, we obtain 

^~0, (v^oo), J 1 = ^y;[2L 3 9 r (2bVi + pVi)] (5-6) 

Therefore, the quantity dJ\/dv tends to be zero at the late time even if we consider initial data 
with Ii ^ 0. So, we cannot show the instability of extreme black holes against electromagnetic 
perturbations by the same way as scalar fields. We may be able to find instability in higher 
order derivatives d™ifi (n > 2). However, since the coupling terms are 0(r 2 ) in the equation, 
we cannot neglect these terms in the limit of r — > when we consider the higher order radial 
derivatives of Eq. (12.16j) . It seems to be difficult problem to show the instability taking into 
account the coupling terms. 



6 Gravitational perturbations 

6.1 Gauge invariant variables on the horizon 

Finally, we study gravitational perturbations. We consider perturbation of GHP variables 
as flij —> Qij + Clij, i&ijk — > ^ijk + ^ijki etc. Here, variables with tildes represent first order 
perturbations. Variables without tildes are background variables. The number of physical 
degrees of freedom of the gravitational perturbations is d(d — 3)/2. On the other hand, 
the number of components of f2y is also d(d — 3)/2. Thus, we can expect that Qij has all 
physical degrees of freedom of the gravitational perturbations. The perturbation variable 0^ 
is transformed by gauge transformations as follows: 
Coordinate transformations (x M — > + £ M (x)): 

n i . q . . + . . . (6.i) 

Spins (tij G so(d — 2)): 

Qij — > Vtij + 2it{i\}Slk\j) ■ (6.2) 

Boosts: 

Qij — > Qij + 2aQ{j . (6-3) 

Null rotations: 

Ci^ -> - 2z k (^ (l S j)k + * (ii)fc ) . (6.4) 

Here, t^, a and z k are infinitesimal functions depending on spacetime coordinates. In 
Eq. (l3.9p . we obtained Q^ = 0(r 3 ) and ^ijk = 0(r 2 ). Thus, under spin and boost transfor- 
mations, we have Q^ — > Q^ + 0(r 3 ). On the other hand, under coordinate transformations 



9 



and null rotations, we have f2jj &ij + 0(r 2 ). Therefore, Q^- [ r=0 and d r £lij\ r= o are gauge 
invariant, but d 2 Clij\ r=0 is not gauge invariant in general. Thus, even if we could show that 
d^Clij\ r= Q blows up along the horizon by the similar way as scalar field perturbations, we 
cannot determine if the instability is physical one or just a gauge mode. We will avoid this 
problem by assuming that the background geometry is algebraically special in sectionJfTH 



6.2 Decoupled equations on the horizon 



We consider the first order perturbation of Eq.f l2.17l) . The right hand side contains cou- 
pling term between f^- and other perturbation variables. We can see that all terms in the 
right hand side are 0(r 2 ). For example, we have (r'p^/)~ = f'p^B + r'p^/ + r'p^ = 0(r 2 ) 
since p and \& are the second order in r. Thus, in Eq. fl2.17p and its radial derivative, the 
right hand side is zero at the horizon. In the left hand side, there is another coupling term, 
4K/ c p / (\I / (jj)fc + ^u8j)k) which can be expanded as 



[4K fc b'(* 



(ij)k 



(6.5) 



This coupling term is 0(r). (Recall that the GHP derivative p' contains the radial derivative 
d r . Hence, we have \)'^ij k = 0(r).) Such a coupling term is harmless when we construct 
a conserved quantity at the horizon and we can show the non-decay of the gravitational 
perturbations by the same way as scalar and electromagnetic perturbations. However, when 
we prove that the perturbations blow up along the horizon, the coupling term is problematic 
since we differentiate Eq. (12.17j) by r. In sect ion J6~4"| we will see that this problem can also be 
avoided by assuming that the background geometry is algebraically special. 



6.3 Conserved quantity on the horizon 



The coupling terms in Eq. (12 . 1 7j) is 0(r) and negligible near the horizon. Thus, near 
horizon, the equation becomes 



d v [2L(2\>% + p%)] = A 2 n ij + 0(r) . 
where the operator Ai is defined as 

•4Aj = --^V fc (L 6 V fc %) +U + Sitfrm - -jy - 2{d - 4)AL 2 J fi, 



(6.6) 



+ 21? f R(i\ k + Rb(i\ k j (l k \j) — 2L 2 R ik jiVL k i 



idk) {l]k - j^(d(L 2 ) A + 2(k - d(L 2 )) (A V k - 2{k - d(L 2 )) k V {l \ 



>.7) 



n 



k\j) 



The derivation of the equation is written in appendixjB] For axisymmetric perturbations 
rrij = 0, the operator A2 is self-adjoint with respect to an inner product 



(Y\Y 2 )= [ 



(6.8) 
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Hereafter, we focus on axisymmetric perturbations. We assume that operator A 2 has a zero 
eigenvalue. In section ]?T2"} we will see that many extreme black holes satisfy this assumption. 
We denote the eigenfunction for the zero eigenvalue as Y^. Operating (Y, *) to both side of 
Eq. (l6.6p . we obtain 

IT p 

= , h = J h Y£ [2L 5 (2b'% + p%)] . (6.9) 

where we used (Y, A 2 Cl) — (A 2 Y, 0) = 0. Therefore, J 2 is a conserved quantity along the 
horizon. Thus, If J 2 7^ at an initial surface, d r Clij and do not both decay along the 
future horizon as v — > 00. Recall that both of d r Clij and fi^ are gauge invariant at the 
horizon. 

6.4 Null rotation to a multiple WAND 

As explained in section JBTT1 and I6T2"} we require conditions ^ijk = 0(r 3 ) and £1^ = 0(r A ) to 
show the instability of the gravitational perturbations. To satisfy these conditions, we assume 
that the background geometry is algebraically special. Then, there is a null rotation (I2.7D 
which transforms the null vector £ to a multiple WAND, that is, ^ijk = Qij = 0. In the near 
horizon geometry (13. 4ft . we have already known the multiple WAND: £ NH = L~ l {2dy + t 2 8r + 
2rk I d I ( j ) ). We can expect that, in the near horizon limit: r = eR, v = V/e and e — > 0, the 
multiple WAND in the full geometry coincides with the £ NH modulo boost transformations^ 
The null vector £ defined in Eq. (13.61) satisfies this condition by itself. (I — > e£ NH in the near 
horizon limit.) Thus, z\ in the null rotation (12.71) should be 0(r 2 ). It follow that the near 
horizon expressions of GHP variables ( 13. 8 j) are correct even after the this null rotation. 

6.5 Instability against gravitational perturbations 

We assume that Cl^ and its tangential derivatives along the horizon decay along the 
horizon. Then, at late time, the conserved quantity I2 becomes 

J 2 ~ 2 f L%*d r Q i3 , (v 00) . (6.10) 

Now, we differentiate Eq. (|2.17j) by r. Since we assumed that the background geometry is 
algebraically special, we can neglect the coupling term in the equation. Near the horizon, the 
equation can be written as 

d v [2L9 r (2b'% + p'Vtij)] = (A 2 + 2)9 r % + V 2 VL i3 + 0(r) , (6.11) 

where, in the linear operator V 2 , there is no radial derivative. Thus, we have T> 2 flij — > 
(v —> 00). The derivation of above equation is in appendixjBl Operating (Y, *) to both side 
of Eq. (l4.6p . we obtain 

-^--h, (v-+oo), ' h = j H Y ii [2^ 5 ^(2b'% + p%)] (6.12) 

2 We checked that this is correct for Kerr-NUT-AdS spacetimes [17j . That is. the multiple WAND in 
Kerr-NUT-AdS spacetimes found in Ref. [T5] approaches £ NH in the near horizon limit. 
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Therefore, the quantity J 2 blow up linearly in time v as 

J 2 ~ hv , (v — > oo) . (6.13) 

Thus, either dfaij or d r flij blows up along the horizon. This implies the instability of the 
extreme black holes against gravitational perturbations. 



7 Unstable extreme black holes 

7.1 Summary of our statement 

Our statements obtained in this paper are as follows: If the operator A s has a zero 
eigenvalue for an axisymmetric perturbation and the horizon conserved quantity I s is non- 
zero, d r ip s and ip s do not both decay along the future horizon as i; — » oo, where ipo = <p, 
ipi = ifi and ip 2 — The explicit expressions of A s are given in Eqs. (l4.2p . (15. 2p and ( 16. 7\i . 
The horizon conserved quantities (14.41) . (15. 4p and (16. 9p are written as 

I s = [ Y*. [2L 2s+1 (2Ws + p'rj>,)] , (7.1) 
Jh 

where Y is the eigenfunction satisfying A S Y = 0. In the proof of this statement, we used 
assumptions ( 13. 3p . 

Hereafter, we assume that ip s and its tangential derivatives along the horizon decay as 
v — > oo. (Then, d r ip s cannot decay.) For scalar field perturbations, we can show that either 
dl4> or d r (p blows up along the horizon. For gravitational perturbations, when the background 
geometry is algebraically special, either d^Qij or d r Qij also blows up along the horizon. For 
electromagnetic perturbations, we could not find the instability in d^cpi or d r <fi. (There may 
be instability in the higher order derivative by r.) 



7.2 Eigenvalues of A s 

The existence of a zero eigenvalue for the horizon operator A s is crucial for the proof of the 
instability. Surprisingly, in study of perturbations of near horizon geometries, the eigenvalues 
of A s have been calculated for some extreme black holes: 4-dimensional extreme Kerr black 
holes [T9 | [20 ] . all 5-dimensional black holes with two rotational symmetries for A = (21] and 
Myers- Perry(-AdS) black holes with equal angular momenta [16 1 122 ) . In this subsection, we 
investigate the existence of zero eigenvalues of A s using their results. 

For massless scalar field, the existence of the zero eigenvalue is trivial since we have 
AqY = when Y is a constant. For massive case, there is no zero eigenvalue in general. In 
some cases, however, it has a zero eigenvalue depending on mass and background geometry. 
For example, in odd-dimensional Myers-Perry- AdS spacetimes with equal angular momenta, 
it was shown that the eigenvalue \ Q is given by 



A Ak(k + N) o , \ / % 

7^ = ^^ L +p\ k = 0,1,2,--- . 7.2 

L A r± 
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where r + is the horizon radius and L 2 is a constant. The integer N relates to the spacetime 
dimension as d — 2N + 3. From this expression, when the scalar field mass is given by 
fi 2 = —4k(k + N) / r^_, the operator Aq has a zero mode. (If the horizon radius r + is sufficiently 
large, the fi 2 does not violate the Breitenlohner-Freedman bound.) 

For electromagnetic and gravitational perturbations (s = 1,2), the existence of a zero 
eigenvalue is not obvious. In 4-dimensional extreme Kerr geometry, the operator A s (s = 
1,2) does not have zero eigenvalue. However, if we consider perturbation equations for Q'- 
instead of flij, we can show the instability [7J. (In the 4-dimensional Kerr geometry, Q'- 
satisfies a decoupled equation although it does not decouple in Myers-Perry spacetimes even 
if we consider the near horizon limit.) For all 5-dimensional black holes with two rotational 
symmetries, eigenvalues \ s are written as^| 

\ 1 = £(£+1) , (£ + l)(£ + 2) , (£ + l)(£ + 2) , 

A 2 = (£ - 1)(£ + 2) , £(£ + 3), £(£ + 3) , (£+ 1)(£ + 3) , (£ + 1)(£ + 3) . 1 ' ' 

where £ = 0, 1,2, We can find that Ai and A2 can be zero for £ = and £ = 0, 1, 
respectively. Thus, in these spacetimes, the gravitational and electromagnetic perturbations 
do not decay in general. In particular, for 5-dimensional Myers-Perry black holes, we can show 
that gravitational perturbations (either d 2 Clij\ r=0 or d r Clij\ r=Q ) blow up along the horizon since 
the spacetimes are known to be algebraically special [23J. 

The horizon induced metrics of Myers-Perry black holes with equal angular momenta can 
be viewed as Hopf fibration over CP N where N is the integer part of (d — 3)/2. Thus, the 
eigenfunction of A s can be decomposed into tensor, vector and scalar harmonics on the base 
space CP N . All eigenvalues for axisymmetric modes are given in Refs. [TB|I2^]. In scalar 
modes, we can always find zero eigenvalues. (Tensor and vector modes can also have zero 
eigenvalues depending on the spacetime dimension d.) Thus, Myers-Perry black holes with 
equal angular momenta in all dimensions are unstable against gravitational perturbations. 
(Electromagnetic perturbations do not decay at least.) Therefore, as far as we calculated, all 
vacuum extreme higher dimensional black holes with vanishing cosmological constant have 
zero eigenvalues in the horizon operator A s . It would be nice if we can show the existence of 
the zero eigenvalues for general black holes. 

8 Discussions 

We studied perturbations in general extreme black hole spacetimes in all dimensions. 
We found a sufficient condition for instability which is summarized in sectionJTJT] Using 
the condition, we showed that 5-dimensional extreme Myers-Perry black holes are unstable 
against gravitational perturbations. For d > 6, we also found gravitational instability in 
extreme Myers-Perry black holes when they have equal angular momenta. 

In the study of perturbations of near horizon geometries [TB|l2ip22] . they considered dimen- 
sional reduction of perturbation equations and obtained effective equations of motion in AdS2- 

3 Note that we need to shift the eigenvalues when we use the results in Refs. [16,21,221 for gravitational 
perturbations since the operator Ai relates to defined by them as C 1 - 2 ' = A2 + 2 for axisymmetric 
perturbations. 
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Then, they used a criterion m 2 < —1/4 to determine the instability of near horizon geometry, 
where m is the effective mass in the AdS2, since this implies violation of the Breitenlohner- 
Freedman (BF) bound. The effective mass relates to eigenvalue of A s as m 2 = \ , Ai, A2 + 2 
for scalar, electromagnetic and gravitational perturbations, respectively. We can see that the 
condition of the instability obtained in this paper (A s = 0) differs from theirs. This discrep- 
ancy comes from difference of type of instabilities. The violation of BF bound (m 2 < —1/4) 
is considered as a condition for an exponential grow of the perturbations. This was explicitly 
shown for scalar field perturbations [T6|. On the other hand, our condition X s = gives a 
power law grow of the perturbations, which is more modest than the exponential one. For 
(d > 6)-dimensional extreme Myers-Perry black holes with equal angular momenta, it was 
shown that the BF bound in the near horizon geometries is violated for gravitational pertur- 
bations |16[l22j. In fact, in the case odd-dimensions, such an instability has been found in 
the full geometry near the extremality [24J. So, the power law instability found in this paper 
may not be important for these spacetimes. However, for 5-dimensional extreme Myers- Perry 
black holes with equal angular momenta, there is no violation of the BF bound in the near 
horizon geometries. In addition to that, from the study of perturbations of full geometries, 
strong evidence of stability for non-extreme black holes has been found in Ref. |25] . Thus, 
the power law instability found in this paper can be important for this spacetime. 

In this paper, we considered only vacuum black holes. Thus, our instability condition 
does not apply to gravitational and electromagnetic perturbations of Reissner-Nordstrom 
(RN) or Kerr-Newman black holes. Further work needs to be done to study their stability 
in extreme limit. They are solutions of N = 2 supergravity and extreme RN black holes 
are supersymmetric. For RN black holes, it is known that the perturbation equations are 
decoupled. Using the decoupled equations, we may be able to find conserved quantities on 
the horizons and show the instability |26j . For Kerr-Newman black holes, the decoupling of the 
perturbation equations has not been succeeded. Hence, even for the non-extreme case, their 
stability has not been studied. As we did in this paper, however, if the coupling terms in the 
perturbation equations are sufficiently small near the horizons, we can study the instability. 
It would be interesting to estimate the order of the coupling terms and study the instability of 
extreme Kerr-Newman black hole. It would make good progress in understanding of stability 
of Kerr-Newman black holes. 

One of the most interesting problems on the instability is its final state. We need to 
solve the time evolution of the instability taking into account the backreaction to specify the 
final state. For scalar field perturbations of RN black holes, we can find an instability for 
spherically symmetric modes [I]. Thus, the evolution equations of the instability are given by 
(1 + 1) -dimensional partial differential equations even if we consider the backreaction. Solving 
the PDEs and finding the final state would be another direction of the future research. 
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HiDGR. 



A GHP variables for extreme black holes 

We take the null basis as in Eq. (13. 6p . The dual one-forms for these vectors are written as 



L 



e 



- ei = —dv , e 1 = e = -r 2 LFdv + 2Ldr , e* = e { = e { - rWdv . (A.l) 



2 



Using the Cartan equations de a + oo ca be c A e b = 0, we can calculate the spin connections u a b c . 
From the definition of L a b in Eq. fl2.2p . we have L a b = —ooboa- Thus, GHP variables for the 
background spacetime are given as 



and 



2r 2r k r 2 F , 

2r 2 Fi L ti {rh a )'^ a 

Ki ~ ~lt~ ' Ti ~~i: + 2L 2 ' 

- — rp , #f' - n t' - -— - ( r/ ^)' g « 

Pij - 2L { j) ' ~~ L 2L 2 



r 2 F)' 2r/i i L* (r/i Q )'e* 



y 10 _ t 1 F? ' 



H r0 , -^n — , Li, 



L 2 ' xx " 2L 2 

2r 2r r 2 F 

Mj = -h a tb e i]a + — ^eg(e i]a ), fc + — ej*(e i]a )' , (A.3) 



where uj ki j is defined by de-i + cukij& k A e J = 0. 



(A.2) 



B Derivation of near horizon equations 

Here, we derive near horizon equations (14. ip . (14. 6p . (15. ip . (15. 5p . ( 16. 6 p and (16. lip . The 
equations for scalar, electromagnetic and gravitational perturbations are written in a unified 
form as 

(2Kb + p'Ws + Bsips = , (B.l) 
where vpo = <p, ipi = tpi and ip2 = ^ij- The angular operators B s are defined by 

£ ^o = (ZA ~ 2t& + pb' - /i 2 )0 , (B.2) 

9/7—3 

S^i = (5,3, - Arjbj + $ - ^3T"A)^ + (-2r^ + 2^ + 2$g + , (B.3) 

2d 

B 2 ip2 = (o fc a fc - 6r fc 3 fc + 4$ - — yA)% 

+ 4(r fc 3 (l - r (i 3 fc + + 4$^)%) + 2§ lkjl tt kl . (B.4) 
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In B s , there is no the radial derivative d r . Up to second order in r, the GHP derivative p can 
be written as 



2 r l 
pV> s = T [d v + r{ikm - b)]ijj s + — d r ip s + r 2 Ctp s + 0(r 3 ) 



(B.5) 



where C is a operator in which the radial derivative d r is not contained. From Eg. ( 13. 131) . we 
obtain [p',r] = 1/(2L). Thus, we have 



(2p'p + p'Ws = d v 



-(2bV s + pV s 



+ ^ |2(zA;m - 6)p' + -j-«9 r + C J V' s + 0(r 2 ) . 



(B.6) 



where C = C + p'{ikm — b). We expand the operator B s as B s = Bf + rB\ + (9(r 2 ). Then, 
Eq. flB.ll) is written as 



Or 



-(2pV s + P> s ) 



+ j^(ikm-b)iff a + Bfiff a 



+ ^ |2(zA;m - 6)p' + + C" J ^ s = ^(r 2 ) , 

where C" = C + LB\/2. Thus, from Eq. (lB~7l) . we obtain 

8 V [2L(2bV s + P> s )] = + 0(r) , 

where 

A, = -2(ikm - b) - L 2 Bf . 
Eq.( lB.8l) expresses Eg. ( 14.in . ( 15. ip and ( 16. 6p . Differentiating Eq.( lB.7p by r, we have 

[2L«9 r (2b> s + p> s )] = [A + 2(6 - 1) - 2ifcm]^« + C '"^ + ( r ) ' 



(B.7) 



(B.8) 
(B.9) 

;b.io) 



where C" = C" + 2(ikm — 6)(p' — d r /(2L)). (Note that there is no radial derivative in the 
operator p' — d r /(2L).) Setting m; = in above equation, we obtain Eqs.( l4.6j ). ( 15. 5 P and 
(16. lip . The explicit expressions of A s can be obtained using near horizon expressions of GHP 
variables and derivatives (13. 8p . (13. lip and ( I3.13p . 



C Useful GHP equations 

We summarize the useful GHP equations for Einstein spacetime satisfying R^ u = Ag^. 
These equations are firstly derived in Ref. [8]. 
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C.l Newman-Penrose equations 

From the Ricci equations, [V M , V^jV^ = R^upaV 7 , we obtain following equations. 

bpij - 6jKi = -pikPkj - Kirf - TiKj - Qij, (NP1) 

br 4 - y*i = Piji-Tj + rj) - (NP2) 
2d\j\Pi\k] = 2np\jk] + 2Kip[ jk] - W ijk , (NP3) 

A 
d~l' 



VPij ~ 9 i r i = ~ T i T j - K i K 'j ~ PikPkj ~ ®ij ~ ~ A 7<%- (NP4) 



Another four equations can be obtained by taking the prime ' of these four. 

C.2 Bianchi equations 

From Bianchi equations, 'V\\O lw \ pa \ = 0, we obtain following equations. 
Boost weight +2: 

b% fe - 23 [j fi fc]i = (2<f>nj\5 k]l - 25u$f k - <$>u jk )Ki 

- 2(*y|(J fl + + + * b -|fl)pi|fc] + 2%t£ ]5 (Bl) 

Boost weight +1: 

-b$y - 8^ + b'% = - (ttjfo - % ik )K k + + 2$f fc + ®5 ik )p kj 

+ - Vi6 jk y k - 2(^ {i 6 j)k + V(ij)k)n - tt ik p' kj i (B2) 
-}>$ ijkl + 28 [fe %,- = - 2*j i | fcI « y] - 2^ m K {l] 

+ ^fjPiki] ~ 2$[ k \iPj\i] + 2$[ k \jpi\i] + 2<& ij[k \mpm\i\ 
- 2^ m T[ 3] - 2%| ij T|' i] - 2VL m p' m + 2tt j[k p' m , (B3) 
-8[3|*i|fcj] =2$y fc |Pi|i] - 2$i\jp k q + $ im [j fc |p m |i] - 2Qiijp' kl] , (B4) 

Boost weight 0: 

tf%jk - 2d\j\<f>i\ k ] =2(%\Sii ~ V]j\u)Pi\k] + (2$i[A]z - 2<J«$/ fc - $^fc)r ; 

+ 2(* i( J b -|, - ^yiO/oJifc] + (B5) 
-28[i$/ fc] =2%p jk] + ^ [y |/9 { | fc ] - 2^ [iP ; fc] - * l{ij \p' m , (B6) 

-3[fc|$tj|H = - %[kl\Pj\m] + *j[/U|Pi|m] - 2^{ k \ i:j p\lm] 

~ ^i[kl\p'j\ m ] + *j[fci|Pi|m] - 2 ^[k\ijP\lm]- (B7) 

Another five equations are obtained by applying the prime operator to above equations. 
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C.3 Maxwell equations 

From Maxwell equations, dF = d * F = 0, we obtain following equations. 

diifii + bF =r-ipi + pijFij - pF - Kitp'i (Ml) 

28^ - bFy =2^^-] + 2Fp m + 2F [i{kPm + (M2) 

2t>Vi + - 3^ =(2p'u1 - - 2i^r,- - 2Fr, + (2p (y) - p<%) (M3) 

8[iFj fc] =(p[iff jk] + ip{iPjk] (M4) 

A further three equations can be obtained by priming above equations. 

C.4 Commutators of derivatives 

The commutation relations for GHP derivatives are given by 

[p, b'] V;,.,.. =(-tj + r'Ad/l] + b ( + + $) T n ... ls 

s 

+ £ (k; /.; - + <r, - T^rj + 2$£J T h ^ A3 , (CI) 

r=l 

[b, 8<] T kl ,„ ks = - (Ki\)' + r^b + Pji?>j)T kl ... ks + 6 (-7^ + /%p^ + T kl ... ks 

s 



r=l 



3i] T fc ,.. fcs = (2p [ii] b' + 24-jb) T fcl ... fcs + 6 (2p ZN p; b1 + 2$g) T fcl ... 



k\...l...k s 



(C3) 



+ ^ (2p Ml |p; b1 + 2p / fcr[i |p ;b1 + + ^— j-%|fc,<%]« J T > 

r=l ^ ' 

The result for [b', 9j] can be obtained from the prime operation of [b, 3j]. 
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